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Abstract —This study presents a method of estimating battery 
cell core and surface temperature using a thermal model coupled 
with electrical impedance measurement, rather than using direct 
surface temperature measurements. This is advantageous over 
previous methods of estimating temperature from impedance, 
which only estimate the average internal temperature. The 
performance of the method is demonstrated experimentally on 
a 2.3 Ah lithium-ion iron phosphate cell fitted with surface 
and core thermocouples for validation. An extended Kalman 
filter, consisting of a reduced order thermal model coupled 
with current, voltage and impedance measurements, is shown to 
accurately predict core and surface temperatures for a current 
excitation profile based on a vehicle drive cycle. A dual extended 
Kalman filter (DEKF) based on the same thermal model and 
impedance measurement input is capable of estimating the 
convection coefficient at the cell surface when the latter is 
unknown. The performance of the DEKF using impedance as 
the measurement input is comparable to an equivalent dual 
Kalman filter using a conventional surface temperature sensor 
as measurement input. 

Index Terms —Lithium-ion battery, impedance, temperature, 
thermal model, Kalman filter, state estimation. 


I. Introduction 

T HE sustainable development of transportation relies on 
the widespread adoption of electric vehicle (EV) and 
hybrid electric vehicle (HEV) technology. Lithium-ion bat¬ 
teries are suitable for these applications due to their high 
specific energy and power density. However, their widespread 
deployment requires reliable on-board battery management 
systems to ensure safe and optimal performance. In partic¬ 
ular, accurate on-board estimation of battery temperature is 
of critical importance. Under typical operating conditions, 
such as a standard vehicle drive cycle, cells may experience 
temperature differences between surface and core of 20 °C 
or more (l). High battery temperatures could trigger thermal 
runaway resulting in fires, venting and electrolyte leakage. 
While such incidents are rare 0 , consequences include costly 
recalls and potential endangerment of human life. 

The conventional approach to temperature estimation is to 
use numerical electrical-thermal models a, a, 0, a, q. 
Such models rely on knowledge of the cell thermal proper¬ 
ties, heat generation rates and thermal boundary conditions. 
Models without online sensor feedback are unlikely to work 
in practice since their temperature predictions may drift from 
the true values due to small uncertainties in measurements and 
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parameters. However, using additional online measurements - 
typically of the cell surface temperature and of the temperature 
of the cooling fluid - coupled with state estimation techniques 
such as Kalman filtering, the cell internal temperature may be 
estimated with high accuracy 0 , a, a, 0 . However, large 
battery packs may contain several thousand cells a, and so 
the requirement for surface temperature sensors on every cell 
represents substantial instrumentation cost. 

An alternative approach to temperature estimation uses 
electrochemical impedance spectroscopy (EIS) measurements 
at one or several frequencies to directly infer the internal cell 
temperature, without using a thermal model a, im ed, 
na, ns. This exploits the fact that impedance is related 
to a type of volume averaged cell temperature, which we 
define later in this article. For brevity, we refer to the use 
of impedance to infer such a volume-averaged temperature as 
‘Impedance-Temperature Detection’ (ITD). This has promise 
for practical application, since methods capable of measuring 
EIS spectra using existing power electronics in a vehicle or 
other application have been developed 01, 03, El- How¬ 
ever, just like conventional surface temperature sensors, ITD 
alone does not provide a unique solution for the temperature 
distribution within the cell. Our previous work showed that 
by combining ITD with surface temperature measurements 
the internal temperature distribution could be estimated ED- 
However, this approach still requires each cell to be fitted 
with a surface temperature sensor. Moreover, whilst the ITD 
technique was validated under constant coolant temperature 
conditions, the accuracy of the technique may be reduced if 
the temperature of the cooling medium is varied rapidly, as 
discussed in Section [IV] Thus, if the cell electrical/thermal 
properties are known or can be identified, this information 
can be exploited to improve the estimate of the thermal state 
of the cell or reduce the number of sensors required. 

In this study we demonstrate that ITD can be used as the 
measurement input to a thermal model in order to estimate 
the cell temperature distribution. First, an extended Kalman 
filter (EKF) is used to estimate the cell temperature distribu¬ 
tion when all the relevant thermal parameters, including the 
convection coefficient, are known. The thermal model consists 
of a polynomial approximation (PA) to the ID cylindrical heat 
equation. The measurement input consists of the cell current 
and voltage, along with periodic measurements of the real part 
of the impedance at a single frequency. Second, a dual ex¬ 
tended Kalman filter (DEKF) is used to identify the convection 
coefficient online when the latter is not known. The predicted 
core and surface temperatures in each case are validated 
against core and surface thermocouple measurements, with 
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agreement to within 0.47 °C (2.4% of the core temperature 
increase). The performance of the combined thermal model 
plus ITD estimator is comparable to the performance of 
the same thermal model coupled with conventional surface 
temperature measurements. Table |T| shows the present study 
in the context of existing temperature estimation techniques, 
which highlights that this is the first study to use impedance 
as the measurement input to a thermal model. 


Study 

Model 

Measurement 



Tsurf 

ITD 

Forgez et al. 0 

/ 

/ 


Lin et al. |[5|, (7| 

/ 

/ 


Kim et al. (6), (4) 

/ 

/ 


Srinivasan et al. 1121, 191 



/ 

Raijmakers et al. 1111 



/ 

Richardson et al. 111711 


/ 

/ 

Present study 

/ 


/ 


Table I 

Comparison of online temperature estimation techniques. 


II. Measurement Principle 

The electrochemical impedance, Z (oj) = Z' ( uj)+jZ" (w), 
of lithium ion cells is a function of temperature, state of 
charge (SOC) and state of health (SOH). Within an appropriate 
frequency range, however, the dependence on SOC and SOH 
is negligible and the impedance can thus be used to infer 
information about the cell temperature lfl2l . Previous ITD 
studies have used as a temperature-dependent parameter the 
real part of the impedance at a specific frequency ITOl , the 
phase shift at a specific frequency Ifl2l . |[9| , and the intercept 
frequency E0- To demonstrate our technique we use the real 
part of the impedance at / = 215 Hz. Our previous work 
showed that the real part of the cell admittance (the inverse of 
the cell impedance) at 215 Hz can be related to the temperature 
distribution using a second order polynomial fit. For an annular 
cell with inner radius ry and outer radius r 0 the real part of 
the admittance is given by 02): 

Y’ = \ f r (ai + a 2 T(r) + a 3 T 2 (r)) dr (1) 
% 7 n 

where ai, a 2 and a 3 are the 1st, 2nd and 3rd coefficients of 
the polynomial relating impedance to uniform cell temperature 
(Y 1 = at + a 2 T uniform + a 3 Tl niform ), provided that the 
admittivity varies in the radial direction only. This assumption 
is valid if the heat transfer from the top and bottom ends of the 
cell is negligible, which is approximately true for cylindrical 
cells connected in series with identical cells on either end 
ED, a configuration which may apply to the majority of 
cells in a large battery pack. However, the application of 
this approach to cooling configurations involving substantial 
end cooling would require a more involved expression for the 
admittance than eq. [I] as well as an appropriate modification 
to the ID thermal model described in the following section. 
Moreover, although the method is applied to a cylindrical cell, 
the proposed approach could be applied to other geometries 
in a similar fashion. The polynomial fit (Fig. [TJ was obtained 



Figure 1. Polynomial fit to experimental data of admittance at / = 215 Hz 
vs. uniform cell temperature. 


by offline impedance measurements on the cell at multiple 
uniform temperatures ED. 

ITD can be viewed as identifying an ElS-based volume 
average temperature Teis , which is defined as the uniform 
cell temperature that would give rise to the measured EISQ 
Thus the impedance input is similar to a conventional temper¬ 
ature measurement since it is a scalar function of the internal 
temperature distribution but does not uniquely identify the 
temperature distribution. Either measurement can therefore be 
used in conjunction with a thermal model, as shown in Fig. 
[2] to estimate core temperature. 
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Figure 2. Schematic of (a) the conventional approach to temperature 
estimation and (b) the proposed approach based on ITD. 


III. Thermal-Impedance Model 
A. Thermal Model 

The cell thermal model consists of the heat equation for ID 
unsteady heat conduction in a cylinder, given by the following 
Boundary Value Problem (BVP) [2]: 

mr.t) , a 2 T(r,t) , k,&T(r,t) , Q(t) 

= + + <2a) 

where p, c p and kt are the density, specific heat capacity and 
thermal conductivity respectively, Vb is the cell volume, and Q 
is the heat generation rate. The boundary conditions are given 

^ote that, since the impedance temperature relationship is non-linear (as 
demonstrated in 1191 ), the ElS-based volume average temperature, 
is not necessarily equal to the volume average temperature, T. Although, it 
should also be noted that the two are typically close in value, particularly if 
the temperature variation within the cell is small. 
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by: 


dT(r,t ) 


dr 
dT(r,t ) 


<9r 


r—0 


— j-(T(r 0 , t) — Too^)) (2b) 
kt 

0 (2c) 


where is the temperature of the heat transfer fluid, and h is 
the convection coefficient. A commonly employed expression 
for the heat source in a lithium ion battery is 


Q = i(v- vocv) + iT°-kssy 


(3) 


which is a simplified version of the expression first proposed 
by Bernardi et al lf20j . The first term is the heat generation 
due to ohmic losses in the cell, charge transfer overpotential 
and mass transfer limitations. The current I and voltage V for 
this expression are measured online. The open circuit voltage 
Uocv is a function of SOC but is approximated here as a 
constant value measured at 50 % SOC, since the HEV drive 
cycles employed in this study operate the cell within a small 
range of SOC (47 — 63 %) and therefore OCV variation. If 
necessary, an estimator of Uocv could also be constructed 
(for example using a dynamic electrical model 1121 1 ). but for 
clarity and brevity we neglect this here. The second term, the 
entropic heat, is neglected in this study because (i) the term 
dU avg /8T is small (0 < dU avg /dT <0.1 mVK 1 ) within 
the operated range of SOC 0, and (ii) the net reversible heat 
would be close to zero when the cell is operating in HEV 
mode. 


B. Polynomial Approximation 

A polynomial approximation (PA) is used to approximate 
the solution of eq. [2a] The approximation was first introduced 
in @ and is described in detail in that article, although the 
essential elements are repeated here for completeness. 

The model assumes a temperature distribution of the form 

T(r,t) = a(t)+ b(t) + 

The two states of the model are the volume averaged temper¬ 
ature T and temperature gradient 7: 

-2 f r ° 2 [ r ° (dT\ 

T = 2 rTdr, 7 = 2 r ( ) dr ^ 

K Jo ri Jo \or J 

The temperature distribution is expressed as a function of T, 
7, and the cell surface temperature, T sur f: 



_ 1 ^~\ r r 

T(r,f) = 4T sur/ -3T -—^7 


— 15r " 

-18T sur/ + 18T+ —^7 


— 45r 

15 T surf - 15T - -^7 


( 6 ) 


Using 2b the surface temperature can be expressed as 
„ 24fc* — 15 ktr„ r n h „ 

surf = 24 k t + r Q h + 48k t + 2 r a h 1+ 24 k t + r 0 h 00 (?) 


By obtaining the volume-average of eq. [2a] and of its 
partial derivative with respect to r, a two-state thermal model 
consisting of two ODEs is obtained: 

x = Ax + Bu 

( 8 ) 

y = Cx + Du 


where x = [T 7 ] , u = [Q T 00 ] T and y = [T core T surf ] T 
are state, inputs and outputs respectively. The system matrices 
A, B, C, and D are defined as: 


A = 

B = 

C = 

D = 


— 48ah —15 ah 

r 0 (24:kt-\-r 0 h) 24kt+r 0 h 

—320a: ft -120a(4fe t +r o fe) 

r%(24k t -\-r 0 h) r^(24k t -\-r 0 h) 

a 4 8cth 

k t Vb r 0 (24k t +r 0 h) 
n 320a/i 

r^(24kt-\-r 0 h) 

24k t -3r 0 h 120r o fc t + 15rg/i 

24kt~\-r 0 h 8(24 fct+r 0 /i) 

24 kt 15 r 0 kt 

24k t -\-r 0 h 48fc t +2r 0 /i 


o 

o 


4 r Q h 

24kt4~r 0 h 

r 0 h 

24 kt-\-T 0 h 


where a = k t /pc p is the cell thermal diffusivity. 


(9) 


C. Impedance Measurement 

Eq. [T] applies to an annulus with inner radius r, and outer 
radius r 0 . If the inner radius is sufficiently small, the cell may 
be treated as a solid cylinder, and eq. |T| becomes 

2 r r ° 

Y' = — r(a 1 +a 2 T(r) + a 3 T 2 (r))dr (10) 

r o Jo 

Substituting eq. [6] in eq. [TO] the real admittance can be 
expressed as a function of of T sur f, T, and 7 

Y' — a\ + a 2 T + 3a 3 T + 2a 3 T 2 ur f — Aa 3 TT sur f 

15a 3 rg7 2 15a 3 r 0 T7 15a 3 r 0 T SMT ./7 00 

+ 32 + 8 8 

Noting from eq. [ 7 ] that T sur f is itself a function of T, 7 and 
Tqo and the cell parameters, admittance is ultimately a function 
of T and 7 . along with the known thermal parameters and 
environmental temperature. In other words, for known values 
of r 0 , k t , c p , p, and h, the impedance is a function of the cell 
state and Too, thus: 

Z' = f(T, 7 , Too) 02) 

IV. Frequency Domain Analysis 

In this section we analyze the error associated with the 
polynomial approximation, by comparing the frequency re¬ 
sponse of the PA model to the frequency response of a 
full analytical solution of eq. ( [2a] ). We also examine the 
approximation employed in 03 , which used a combination 
of impedance and surface temperature measurements but no 
thermal model. To achieve a unique solution in that case, it was 
necessary to impose the assumption of a quadratic solution to 
the temperature distribution, and so we refer to this here as 
the ‘quadratic assumption’ (QA) solution. 
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As in ( 6 ]|, the frequency response function of the PA thermal 
system, H(s), is calculated by 

H(s) = D + C(sl — A) -1 B (13) 

where s = joj and I is the identity matrix. The frequency 
response of the analytical model is derived in l 22 l . and that 
of the QA solution is derived in Appendix [A] 



20 



Frequency (Hz) 



Frequency (Hz) 


Figure 3. Comparison of frequency responses of (i) analytical solution to 
the ID cylindrical heat transfer problem, (ii) the polynomial approximation 
used in the current study and (iii) the quadratic assumption used in ED- 

Fig. [3] shows the impact of changes in heat generation on 
T core and T sur f (Hu and H 21 respectively), and the impact of 
changes in cooling fluid temperature on T core and T sur f 
(H 12 and H 22 respectively), for each model along with the 
error relative to the analytical solution. The results in these 
plots were obtained using the thermal parameters of the 26650 
cell used in the present study (Table 0 - The response of 
Hu and H 21 for both the PA and QA models are in good 
agreement with the analytical solution (note that the error 
in the response of H 21 and H 22 for the QA model is 0 
since it takes the measured surface temperature as one of 
its inputs). However, the responses for both cases to changes 
in Too are less satisfactory. In particular, the response H 12 
for the QA model shows a rapid increase in error relative to 
the analytical solution at frequencies above ~ 1 0“ 3 Hz. The 
PA performs well up to higher frequencies, although its error 
in H 21 and H 2 2 become unsatisfactory above ~ 10 ” 2 Hz. 
However, the frequency range at which the PA agrees with 
the analytical solution is broader than that of the QA model, 
and is considered satisfactory given the slow rate at which the 
cooling fluid temperature fluctuates in a typical battery thermal 
management system. 


V. Experimental 

Experiments were carried out with a 2.3 Ah cylindrical cell 
(A123 Model ANR26650 ml-A, length 65 mm, diameter 26 


mm) with LiFePCH positive electrode and graphite negative 
electrode. The cell was fitted with two thermocouples, one 
on the surface and another inserted into the core via a hole 
which was drilled in the positive electrode end (Fig. [4]). Cell 
cycling and impedance measurements were carried out using a 
Biologic HCP-1005 potentiostat/booster. The impedance was 
measured using Galvanostatic Impedance Spectroscopy with a 
200 mA peak-to-peak perturbation current. The environmental 
temperature was controlled with a Votsch VT4002 thermal 
chamber. The chamber includes a fan which operates contin¬ 
uously at a fixed speed during operation. 

In order to calibrate the impedance against temperature, 
EIS measurements were first conducted on the cell in thermal 
equilibrium at a range of temperatures. These experiments and 
the subsequent identification of the polynomial coefficients a±, 
ci 2 , and (-(3 are described in 021 . 

Dynamic experiments were then conducted using two 3500 s 
current excitation profiles - the first to parameterise the thermal 
model, and the second to validate the identified parameters 
and to demonstrate the temperature estimation technique. The 
profiles were generated by looping over different portions of 
an Artemis HEV drive cycle. The applied currents were in the 
range —23 A to +30 A. For the duration of the experiments, 
single frequency (215 Hz) impedance measurements were 
carried out every 24 s and the surface and internal temperatures 
were also monitored. In order to minimise heat loss through 
the cell ends, these were insulated using Styrofoam (Fig. [4|. 
Before each experiment, the SOC was adjusted to 50% by 
drawing a 0.9 C current. The temperature of the thermal 
chamber was set to 8 °C and the cell was allowed to rest 
until its temperature equilibrated before experiments began. 

Besides being a function of temperature and SOC, the 
impedance is also a function of DC current, mainly due to 
the charge transfer polarization decreasing with increasing 
current ll23ll . Previous results confirmed that when the EIS 
perturbation current is superposed over an applied DC current, 
the impedance measurement is altered ifTTl . To overcome this, 
the cell was allowed to rest briefly for 4 s before each EIS 
measurement was taken, i.e. 20 s periods of the excitation 
current followed by 4 s rests. The duration of this rest period 
was kept as short as possible to ensure that the thermal 
response of the cell to the applied cycle was not significantly 
altered, and it was found that the core cell temperature dropped 
by at most 0.25°C during these rest periods. The issue of 
impedance measurement under DC current warrants further 
investigation. 

VI. Model Parameterisation & Validation 

Parameterisation is performed to estimate the values of 
k t , c p , and h. The density p was identified in advance by 
measuring the cell mass and dividing by its volume. The 
values from the first excitation profile (comprising cell current, 
voltage plus surface, core and chamber temperatures) were 
used for the estimation. 

The parameterisation was carried out offline using finiri¬ 
se arch in Matlab to minimise the magnitude of the Euclidean 
distance between the measured and estimated core and surface 
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Temperature sensors 
(for validation): T core , T surf 



d 



Figure 4. Experimental setup for parameter identification and validation, (a) 
schematic diagram with cutaway view showing cell core and jelly roll, (b) 
cell drilling procedure, (c) prepared cell inside thermal chamber, (d) power 
supply & thermal chamber. 


temperatures, as in a. Table [II] compares the identified 
parameters with the initial guesses and parameters for the 
same cell from the literature. The estimated parameters are 
close to those reported in the literature. The deviations may 
be attributed to manufacturing variability, error in the heat gen¬ 
eration calculation (due to the omission of entropic heating in 


all of these studies), heat generation in the contact resistances 
between the cell and connecting wires and/or measurement 
uncertainty in the temperature. The convection coefficient is 
within the range expected of forced convection air cooling 

EH- 


Parameter 

Units 

Reference 

Initial 

Identified 

p 

kg m'- 3 

2047-2118 

0 ED, CD 

- 

2107 

c p 

J kg- 1 K 1 

1004.9-1109.2 

0, 0. 0 

1050 

1171.6 

kt 

W m 1 K 1 

0.488-0.690 

(22J, EH. 0 

0.55 

0.404 

h 

W m 2 

- 

20 

39.3 


Table H 

Comparison of reference & estimated parameters 


The measured core and surface temperatures (subscript 
‘exp’) and the corresponding model predictions (subscript 
‘m’) for the parameterised model are shown in Fig. [5] The 
root-mean-square errors (RMSE) in the surface and core 
temperatures are 0.19 °C and 0.18 °C respectively. 



Figure 5. Model parameterisation: Comparison between measured and 
predicted core and surface temperatures in the parameterised model. 

The model with identified parameters was validated against 
the second current excitation profile (Fig. [6}. The RMSEs 
in the core and surface temperatures were 0.21 °C and 0.16 
°C respectively in this case. These errors are only marginally 
greater than those in the parameterisation test, indicating that 
the estimation is satisfactory. 



Figure 6. Model validation. Comparison between measured and predicted 
core and surface temperatures in the parameterised model applied to the 
second current excitation profile. 

VII. State Estimation 

Kim et al. 0 showed that the effect of changes to the value 
of h on the predicted surface and core temperatures is greater 
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than the effect of changes to the other thermal parameters. 
Moreover, h depends strongly on the thermal management 
system settings and its calculation often relies on empirical 
correlations between coolant flow rates and heat transfer. 
Thus, there is a need to identify the convection coefficient 
online during operation. This section outlines the use of a 
dual extended Kalman filter (DEKF) ||27| for estimating the 
core and surface temperatures and the convection coefficient. 
The DEKF reduces to an EKF if the convection coefficient is 
assumed known and provided to the model in advance. 

We firstly modify eq. [9] by rewriting it as a discrete 
time model, setting the impedance as the model output, and 
explicitly including the dependency on the parameter, h k : 

x fc+1 = A (/ifc)x fc + B(h k )u k + v fc (14) 

Vk = /(xfc, h k ) +n k (15) 

hk-\-i — h k T e k (16) 

where y k = Z' and f(x. k , h k ) is the non-linear function 
relating the state vector to the measurement (i.e. eq. [13- 
and v^., n k and e k are the noise inputs of the covariance 
matrices R v , R" and R e . The states, inputs and measured 
outputs are thus x = [T 7] , u = [Q T^] and y = Z'. 
Note that, although the impedance is now the model output, 
the core and surface temperatures are also computed from the 
identified states and parameter at each time step using eq. [8] 
for validation against the thermocouple measurements. A and 
B are system matrices in the discrete-time domain, given by 

A = e ( AAt \ B = A -1 (A — I)B (17) 

where At is the sampling time of 1 s. The update processes 
are then given as follows. The time update processes for the 
parameter filter are: 


where K* is the Kalman gain for the state, and H* is the 
Jacobian matrix of partial derivatives of / with respect to x: 


H 


X 

k 


df{x k ,h k ) 

<9xfc 


x fc =x fc 


(26) 


The measurement update processes for the parameter filter are: 
K h k = (P*)"(I#) T (H k (P k )~ (H k ) T + iT)’ 1 (27) 

h k = h k + K k [z k - /(xfc, h~)^j (28) 

P£ = (I - K h k H%){Pj;)- (29) 


where H k is the Jacobian matrix of partial derivatives of / 
with respect to h: 


H 


h 

k 


<9/(xfc,/ifc) 


dh k 


hu—h. 


(30) 


The above algorithm can be simplified to a standard EKF by 
omitting the parameter update processes (eqs. p~6l p~8|n~9| and 
27p0| and assuming the convection coefficient is fixed. In the 
following section we investigate the performance of both the 
baseline EKF and the full DEKF algorithm. 


VIII. Results 

We first investigate the performance of an EKF estimator 
whereby the convection coefficient is provided and assumed 
fixed. We then compare the performance of the DEKF algo¬ 
rithm with that of the baseline EKF when an incorrect initial 
estimate of the convection coefficient is provided. Lastly, we 
compare the performance of the DEKF with that of a dual 
Kalman filter (DKF) based on the same thermal model but 
with T sur f as the measurement input rather than Z'. 


h k = h k _ 1 (18) 

(P' 1 )- = Pti + R e (19) 

where h k and h k are the a priori and a posteriori estimates of 
the parameter h, and {P k )~ and P k _ 1 are the corresponding 
error covariances. 

The time update processes for the state filter are: 

Xfc = Afc_iXfc_i + Bfc_iUfc_i (20) 

(Pfc)“ = A fc _ 1 Pfc_ 1 Afc_i + R v (21) 

where x^ and Xfc are the a priori and a posteriori estimates 
of the state, and (P*) - and are the corresponding error 

covariances. The matrices Afc_i and Bfc_i are calculated by: 

Afc_i = A(h)\ h=h ., Bfc_! = B {h)\ h=h . (22) 

Since the relationship between impedance and the cell state 
is non-linear, the measurement model must be linearised 
about the predicted observation at each measurement. The 
measurement update equations for the state filter are: 

Kfc = (P£)-(H£) T (H£(P£)-(H£f + iT)- 1 (23) 

Xfc = Xfc + [z k - /(Xfc , hi )) (24) 

Pfc = (I — KfcHfc)(Pfc)“ (25) 


A. Convection Coefficient Known 

The initial state estimate provided to the battery is xo = 
[25 0], i.e. the battery has a uniform temperature distribution 
at 25 °C. The true initial battery state is a uniform temperature 
distribution at 8 °C. The covariance matrices are calculated as 
R n = and R v = j3^diag( 2, 2). The first tuning parameter 
is chosen as o n = 1 x 10 -4 fl, which is a rough estimate 
of the standard deviation of the impedance measurement. The 
second tuning parameter was chosen as (3 V = 0.1, by trial and 
error. Fig. [7] shows that, using the EKF, the core and surface 
temperatures quickly converge to the correct values and are 
accurately estimated throughout the entire excitation profile. 
The RMSEs of core and surface temperature are 1.35 °C and 
1.34 °C respectively. In contrast, the RMSEs for the open 
loop model (subscript ‘m’) with no measurement feedback 
are 6.66 °C and 4.42 °C respectively. It should be noted that 
since the uncertainty of the impedance measurement typically 
increases as impedance decreases, the temperature estimates 
could be more uncertain at higher temperatures. Hence, the 
implementation of this technique could be more challenging 
at higher ambient temperature conditions than those studied 
here. 

It should be noted that we also achieved similar performance 
using a simpler EKF based on Z' with the assumption that 
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Figure 7. Temperature results for EKF using Z' as measurement input 


the impedance is related directly to T rather than to Teis- 
However, this assumption may lead to unsatisfactory results for 
cells with a larger radius or when larger temperature gradients 
exist within the cell. Moreover, since this approach assumes 
that the impedance is a function of the state only (and not 
the parameter h ), it is not suitable for the application of the 
DEKF discussed in the following section. 



t/s 


Figure 8. Temperature results for DEKF using Z' as measurement input 


B. Convection Coefficient Unknown 

Next we investigate the performance of the DEKF. The 
same incorrect initial state estimate is provided to the battery, 
xq = [25 0]. Moreover, an incorrect initial estimate for the 
convection coefficient is provided, ho = 2 x h trU e■ This 
value of h is also provided to the EKF. The error covariance 
matrix for the parameter estimator is R e = where the 
tuning parameter is chosen as /3 e = 2.5. Fig. [8] compares 
the results of both of these cases against the thermocouple 
measurements. The EKF is shown to overestimate the core 
temperature and underestimate the surface temperature for 
the duration of the experiments. This is expected, since the 
impedance measurement ensures the accuracy of the volume 
averaged cell temperature but the model assumes that the 
convection coefficient is higher than in reality, and therefore 
the temperature difference across the cell is overestimated. In 
contrast, the DEKF corrects the convection coefficient, and 
thus improves the accuracy of the subsequent temperature pre¬ 
dictions. This is evident from the errors in the core temperature 
estimate (top plot of Fig. [8]), which initially are similar in both 
cases but drop to much smaller values for the DEKF once the 
correct convection coefficient is identified. The RMSEs of core 


and surface temperature in each case are shown in Table III 


along with the values for the time period, 1200 < t < 3500 s 
(i.e. after the convection coefficient value has converged). 

Finally, we investigate the case where T aur f is used as the 
measurement (y = T sur f) to the estimator rather than Z'. 
This results in a linear KF and DKF exactly equivalent to 
that studied in 0. The tuning parameters for the covariance 
matrices are also chosen to be the same as those employed 
in 0. The same initial state and parameter estimates are 
provided as for the DEKF. Fig. [9] shows that the standard 


EKF in this case overpredicts the core temperature by a 
greater margin than the EKF based on Z ', although the surface 
temperature estimate is much more accurate. This is because 
the surface thermocouple ensures an accurate surface temper¬ 
ature estimate and to reconcile this with the overestimated 
convection coefficient, the core temperature estimate is forced 
to be much greater. The DKF correctly identifies the correct 
convection coefficient, in the same way as the DEKF. Since 
the thermocouple measurement exhibits less noise than the 
impedance measurement, the model converges to the correct 
estimate for h more quickly than in the DEFK case, as shown 


by the RMSE values in Table III 


Method 

0 < t < 3500 s 

1200 < t < 3500 s 


Tcore 

Tsurf 

Tcore 

Tsurf 

EKF + Z' 

2.04 

2.06 

1.79 

1.98 

KF + T aur f 

2.49 

1.44 

2.90 

1.58 

DEKF + Z' 

1.43 

1.24 

0.47 

0.42 

DKF + T sur f 

0.36 

0.33 

0.16 

0.14 


Table III 

Comparison of RMSEs for core and surface temperatures (°C) 

WITH UNKNOWN CONVECTION COEFFICIENT. 


In conclusion, the temperature and convection coefficient 
estimators using Z' as measurement input are capable of ac¬ 
curately estimating the core and surface temperatures and the 
convection coefficient. The performance is comparable to that 
of an estimator using the same thermal model coupled with 
surface temperature measurements. Moreover, the performance 
of the present method is superior to that of previous methods 
based on impedance measurements alone, which only provide 
an estimate of the average internal temperature of the cell. 
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Figure 9. DKF using T sur f as measurement input. 


IX. Conclusions 


Impedance temperature detection enables both core and sur¬ 
face temperature estimation without using temperature sensors. 
In this study, the use of ITD as measurement input to a cell 
thermal model is demonstrated for the first time. 

Previously, we estimated cell temperature distribution by 
combining ITD with a surface temperature measurement and 
imposing a quadratic assumption on the radial temperature 
profile. Frequency domain analysis shows that the QA solution 
may be inaccurate if the temperature of the cooling fluid has 
fluctuations on the order of 10 -3 Hz or higher. The PA thermal 
model used in the present study is robust to higher frequency 
fluctuations (~ 10 -2 Hz). 

An EKF using a parameterised PA thermal model with ITD 
measurement input is shown to accurately predict core and 
surface temperatures for a current excitation profile based on 
an Artemis HEV drive cycle. A DEKF based on the same 
thermal model and measurement input is capable of accurately 
identifying the convection coefficient when the latter is not 
provided to the model in advance. The performance of the 
DEKF using impedance as measurement input is comparable 
to an equivalent DKF estimator using surface temperature as 
measurement input, although the latter is slightly superior due 
to the higher accuracy of the thermocouple. 

Future work will investigate the application of ITD to 
multiple cells in a battery pack, as well as investigating 
methods of combining impedance with conventional sensors to 
enable more robust temperature monitoring and fault detection, 
and self-calibration. 
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Appendix 


A. Frequency Domain Analysis of Quadratic Assumption 

The analysis leading to the frequency response plots of the 
QA model in Fig. [3] is outlined in this section. The QA model 
was used in 113 to obtain a unique solution for the temperature 
distribution when the impedance and surface temperature were 
measured but the cell thermal properties and heat generation 
rates were assumed unknown. To achieve a unique solution 
in that case, it was necessary to impose the assumption of a 
quadratic temperature distribution based on the solution of the 
ID heat equation at steady state. 

Muratori et al. l22l showed that the time domain PDE of 


eq. 2a can be transformed into an equivalent ODE problem in 


the frequency domain which can be solved analytically. Using 
this approach the solution for the temperatures at the core and 
surface of the cell are shown to be: 


1 


T core (s) — - — ^Q(s) + 
k t a z 

Tsurf (s) — - nQ(s) + 

k t a z 


k t {Too (s) fe 3 2 Q(s)) 

Tj 0 (ar o ) - aJi(ar 0 ) 

j^C?oo(s) — -jr ^2 Q(s)) 
Tj 0 (ar o ) - aJ\{ar 0 ) 


Jo(c 


(31) 


(32) 


where a 2 = sa 1 , Q(s) is the transform of Q(t), and Ji 
is the i th - order Bessel function of the first kind (28). Eqs. 


|3T| and 32 can be interpreted as the outputs of a continuous 
time dynamic system (22l , where u(t) = [Q(t), T x {t)] T and 
y (t) = [T core (t), T sur f{t)] T , such that the solution of the 
BVP is equivalent to the impulse response of the system: 


(33) 


where the H matrix is formed by the transfer functions: 


T-'core (^) 


'Hu(j) 

H 12 (s)" 

' Q{s )' 

Tsurf (^)_ 


H 21 (s) 

H 22 (s)_ 

Too(s)_ 


Hu(s) 

H 12 («) 

H 21 (s) 

H 22 (s) 


1 ^J 0 (ar o ) - aJi(ar 0 ) - T 

ha 2 Y t Jo( ar o ) - aJi(ar 0 ) 
h 

_ kt _ 

Y t Jo{ar 0 ) - aJi(ar 0 ) 

1 —aJi(ar 0 ) 

ha 2 Y t J o( ar °) ~ aJi(ar 0 ) 

T t M ar °) 

Tj 0 (ar o ) - aJi(ar 0 ) 


(34) 

(35) 

(36) 

(37) 


This system of transfer functions gives frequency responses 
for the analytical solution results plotted in Fig. [T] Using a 
similar approach, we can develop an analytical solution to the 
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QA model used in HZ), where we denote the new transfer 
function Hq A (s'). In this case, both the volume-averaged cell 
temperature (identified via the impcdancf]) and the surface 
temperature were measured directly, and it is assumed that 
no other information was available. Since these two inputs 
alone are not sufficient to achieve a unique solution for the 
temperature distribution, it was also necessary to impose the 
following constraint on the temperature profile: 

Tqa (f ) = T sur f + ( T corCt q A — T sur f) ^1 - — ^ (38) 


where the H matrix is formed by the functions: 


Hqa, n (s) 


kta 2 r 0 Ji(ar 0 ) — 2hr 0 aJo(ar 0 ) + AhJ\(ar 0 ) 
k t r 0 a 3 ( hJ 0 (ar o ) - fc t aJi(ar 0 )) 


— hk t a 3 Joar 0 + 4hk t a 2 J\ar 0 
QA ' 21 k t r 0 a 3 ( hJ 0 ar o - k t aJi{ar 0 )) 

H QA, i2(s) = H 21 (s) 

HqA,22(s) = H-22 (s) 


(44) 

(45) 

(46) 

(47) 


References 


This is the ID steady-state solution of the heat equation for 
a cylinder with uniform heat generation l24l . with T(r) and 
T core replaced by T QA {r) and T cor ^Q A - 

The volume average of the QA temperature distribution is 
set equal to the true volume average temperature, giving: 

Tq A = \ I* rT QA (r)dr = T (39) 

r o Jo 

Substituting eq. [38] in eq. [39] and integrating, the QA core 
temperature becomes 


Tcore, QA = 2T — T sur f 


(40) 


Thus, an estimate for the core temperature is obtained di¬ 
rectly from the surface and volume averaged temperature 
measurements. Since the surface temperature is measured, 
the values of Hq,a, 2i(s) and Hqa, 22 (s) are identical to 
the corresponding values of the unsteady thermal model. 
The values of Hq^i, u(s) and Hq A . 12 (s) can be obtained 
as follows: Substituting T from eq. [ 5 ] into eq. 40 the QA 
approximation of the core temperature can be expressed as a 
function of the temperature distribution: 


Tcore, QA = \ f (rT(r)dr - T surf ) dr 
r o Jo 


(41) 


Substituting eqs. 31 and [32] for T(r) and T sur f, respectively, 
and integrating (noting that ff r J 0 (ar)dr = r a Ji(ar 0 )/a), 
we obtain: 


T, 


core , QA 


Q(s) 

k t a 2 


^(Too(a)-^Q(a)) 


AJi(ar 0 ) 


T-J 0 {ar o ) - aJi(ar 0 ) 

Thus, the QA system model is given by: 


- J 0 (ar o ) 


(42) 


Tcore, Qa(s) 


Hqa,h(s) 

Hqa,12(s) 

' Q(a) ' 

T S urf, QA (^)_ 


Hqa,2i(«) 

Hqa, 22 (s)_ 

Too{s)_ 


(43) 

2 As discussed in Section firthe impedance actually identifies T eis-> which 
is not necessarily equal to T. However, the assumption that Teis = 
is satisfactory for the purpose of identifying the frequencies at which errors 
relative to the analytical solution become significant, and is only used for this 
purpose in this article. 
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